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It is well established that an entanglement encoded in the Bell states of a two-qubit system with 
correlated spins exhibits completely different evolution properties than that encoded in states with 
the anti-correlated spins. A complete and abrupt loss of the entanglement, called the entanglement 
sudden death, can be found to occur for the spin correlated states, but the entanglement evolves 
without any discontinuity or decays asymptotically for the spin anti-correlated states. We consider 
the evolution of an initial entanglement encoded in the spin anti-correlated states and demonstrate 
that the asymptotic behavior predicted before occurs only in the weak coupling limit or equiva- 
lently when the rotating-wave approximation (RWA) is made on the interaction Hamiltonian of the 
qubits with the field. If we do not restrict ourselves to the RWA, we find that the entanglement 
undergoes a discontinuity, the sudden death phenomenon. We illustrate this behavior by employing 
an efficient scheme for entanglement evolution between two cold-trapped atoms located inside a 
single-mode cavity. Although only a single excitation is initially present in the system, we find that 
the two-photon excited state, which plays the key role for the discontinuity in the behavior of the 
entanglement, gains a population over a short time of the evolution. When the RWA is made on the 
interaction, the two-photon excited state remains unpopulated for all times and the discontinuity is 
absent. We attribute this phenomenon to the principle of complementarity between the evolution 
time and energy, and the presence of the counter-rotating terms in the interaction Hamiltonian. 

PACS numbers: 03.65.Ud, 03.67.-a, 42.50.Pq 



An understanding of entanglement evolution and en- 
tanglement transfer between qubits is of fundamental in- 
terest in quantum information processing [l], [2|. The 
controlled transfer that preserves initial entanglement is 
crucially important. Transfer processes are susceptible 
to decoherence and dissipation due to the inevitable cou- 
pling of the qubits and the transfer channels to an ex- 
ternal environment. Therefore, in order to minimize de- 
coherence effects and to achieve the perfect fidelity, fast 
transfer processes or transfer operations performed over 
a very short time scale are highly desirable. 

It has been recognized that the entanglement evolution 
depends on the state in which it is encoded. For a simple 
system of two qubits, the basis states for entanglement 
are four mutually orthogonal Bell states Q. The states 
can be divided into two groups, one involving linear su- 
perpositions of the spin correlated states and the other 
involving spin anti-correlated states. The states belong- 
ing to these groups are often called two-photon and one- 
photon entangled states, respectively. The qubits can be 
prepared in a spin correlated state by a transfer of two- 
photon entangled states from quantum-correlated light 
fields produced e.g., in a nonlinear process of parametric 
down conversion [4l| . Preparation of a spin anti-correlated 
state is more sophisticated as it involves a single excita- 
tion "shared" between two qubits. In principle, it can 
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be achieved, for example by applying a short single laser 
pulse either in a running or in a standing wave configura- 
tion 0. This will result in the qubits prepared either in 
a symmetric or in an antisymmetric combination of the 
spin anti-correlated states. 

Dynamics of an entanglement encoded in spin corre- 
lated states have been extensively studied since the pio- 
neering work of Yu and Eberly Q , who showed that 
an initial entanglement encoded in two separate qubits 
interacting with local environments can decay to zero in 
a finite time When the qubits are subjected of the 
interaction with each other through the coupling to the 
same environment, the already dead entanglement may 
revival after a finite time Q. The interaction between 
the qubits induces a population difference between the 
symmetric and antisymmetric combinations of the spin 
anti-correlated states, which results in a nonzero entan- 
glement. 

A completely different conclusion applies to the evo- 
lution of entanglement initially encoded in a s pin anti- 
correlated state. It was pointed out by Jamroz flO'l that 
an initial entanglement encoded in a spin anti-correlared 
state of two independent qubits interacting with local en- 
vironments decays asymptotically in time without any 
discontinuity. This prediction agrees with Yonac and 
Eberly 's work [ll[, and also with other results [l^l- The 
same conclusion applies to the case of two qubits mutu- 
ally interacting through the coupling to a common en- 
vironment [Tsj . In these papers, the analysis were re- 
stricted to the RWA and this raises the question on the 
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validity of the results in the strong coupling regime where 
the breakdown in the RWA occurs |14| . 

In contrast to the numerous publications concerning 
atomic dynamics under the RWA, there are only a limited 
number of studies beyond the RWA. In terms of appli- 
cations, the non-RWA calculations exploit short time or 
high-intensity dynamics where interesting effects appear 
not observed when the RWA is introduced. In this con- 
nection, we should mention the work on the Bloch-Siegert 
shift [3], chaos in the Jaynes-Cummings model bi- 
furcations in the phase space [l7| and a fine structure in 
the optical Stern-Gerlach effect In the connection to 
entanglement, correlations between two separate atomic 
ensembles have recently been analyzed by Ng and Bur- 
nett fl^. It has been demonstrated that the ensembles 
interacting with the cavity field and initially prepared in 
their ground states undergo time evolution and become 
entangled over a short time only if the RWA is not made 
on the interaction Hamiltonian. 

In this paper, we consider dynamics of the spin anti- 
correlated states of two identical qubits without the pres- 
ence of any external excitations. We treat the problem 
fully quantum mechanically and do not restrict ourselves 
to the RWA. We find significant quantitative deviations 
from the RWA for the time evolution of an initial entan- 
glement encoded in the spin anti-correlated states. We 
show that the entanglement may undergo a discontinu- 
ity that, on the other hand, requires either an initial or 
a transient buildup of the population in the two-photon 
state of the system. One could argue that this rules out 
the discontinuity in the entanglement evolution since only 
a single excitation was present initially. We shall demon- 
strate that this is not the case if one considers an evolu- 
tion with the non-RWA Hamiltonian and interpret this 
result is a consequence of the principle of complementar- 
ity between the evolution time and energy. 

To describe the entanglement evolution in a two qubit 
system, we use concurrence, an entanglement measure 
that relates entangled properties to the coherence prop- 
erties of the qubits [20|. In order to compute the con- 
currence, one needs the density matrix of the two qubit 
system written in the basis of the product states |1) = 

1 ii), |2) = I Ti), |3) = I it), |4) = I TT). Here, | Ti) repre- 
sents the qubit 1 in the excited ("up") state, the qubit 

2 in the ground ("down") state. The density matrix is 
in general composed of sixteen nonzero elements. We 
make an assumption that the qubits evolve without the 
presence of any external fields. In this case, the density 
matrix takes a simple block diagonal form 



For a system described by the density matrix ([1]), the 
concurrence has the form 



Piiit) 

P22{t) P23{t) 

P32(t) PZ3{t) 

pii{t) 



(1) 



in which we keep all the diagonal elements (popula- 
tions) and the coherences P23{t) and P32{t) that might 
be nonzero initially or can build up during the evolution 
of the system. 



with 



C{t) = max{0, C{t)}. 



C{t) = 2|p23(i)| - 2 v/pii(i)p44(i), 



(2) 



(3) 



from which we see immediately that the discontinuity 
behavior cannot be achieved in any state of the system 
if only a single excitation was present since pAiit) = 0. 
In this case, the concurrence depends solely on the co- 
herence P2'i{t)- The discontinuity or threshold behavior 
of the concurrence requires a nonzero population of the 
ground and the upper two-photon states of the system. 
If an initial excitation were redistributed during the evo- 
lution among all of the states, we would have a possi- 
bility for a discontinuous evolution of entanglement. We 
then would have a prototype of the entanglement sudden 
death. 

To gain further insight into the discontinuity behavior 
of entanglement, we employ the Bell states that are max- 
imally entangled states of two qubits with correlated and 
anti-correlated spins 

k> - ^(in> + uT», ia> = i=(in>-uT», 

H = ^(|TT) + IU», |/3)^-^(|TT)-|U»,(4) 

In terms of the Bell states, the concurrence ([3]) has the 
form 

C{t) - \j [Pss (t) - Paa it)f- [psa (t) - Pas {t)f 



- \J [Paa [t) + Ppp {t)f~ [Pafi {t)+Pl3a {t)f ■ (5) 

As the entanglement has been defined by the require- 
ment that C{t) > 0, it follows immediately from Eq. ([5]) 
that the entanglement depends on the distribution of 
the population between the spin correlated and anti- 
correlated states. An entanglement creation involving 
the spin anti-correlated states is diminished by the pres- 
ence of population and coherence between the spin corre- 
lated states. The evident competition between the spin 
correlated and anti-correlated states in the creation of en- 
tanglement may lead to a discontinuity in the time evolu- 
tion of the entanglement. Therefore it immediately raises 
the question of whether a discontinuity in the evolution 
of the entanglement can ever be achieved if initially only 
a single excitation was present and no external excita- 
tions are applied to the system. This problem is treated 
quantitatively below. 

Let us now examine the above general results for the 
properties of the concurrence on a simple example of two 
identical two-level atoms (qubits) located inside a high-Q 
single- mode standing- wave cavity [2lj. We assume that 
the qubits undergo a strong coupling to the cavity mode 



3 



when located at the antinode of the cavity field, and con- 
sider the situation where the strength of the coupling is 
controlled by varying the position of the qubits inside the 
standing-wave cavity mode. The dynamics of the system 
are determined by the master equation of the density op- 
erator ps of the total, qubits plus the cavity field system 

= -j[H,ps\ - -K {a^aps + PsO^a - 2apsa^) , (6) 
where k is the damping rate of the cavity mode, and 



1 ^ 

i^nRWA = ^ftt^o ^ cTj + fiuia^a 



+nY^[g{r,)a^^a^ +g*{r,)aa^] (7) 

is the non-RWA Hamiltonian of the qubits of transition 
frequency interacting with the single-mode cavity field 
of frequency lo. The operators a and are the annihila- 
tion and creation operators of the cavity field, while <Tx 
and (Tz are the Pauli matrices in the x and z directions, 
respectively. The parameter g{rj) is the coupling con- 
stant between the jth qubit and the cavity mode at a 
position rj of the qubit along the cavity axis. The cou- 
pling constant varies with the position of the qubits as 



9{ri,2) =5osin[7r(iTd)/A], 



(8) 



where L is the size of the cavity and A is the cavity wave- 
length. We assume that the qubits are placed symmet- 
rically about the antinode of the cavity mode such that 
ri + r2 = L and r2 — ri = d is the distance between the 
qubits. When d = 0, the qubits are then at the antin- 
ode of the cavity mode and experience the peak coupling 
strength g^. As d increases, the coupling strength de- 
creases and approaches zero at d = L. 

When the RWA is made, the Hamiltonian ([7]) takes the 
form 



I 2 

■ffRWA = -hujQ^aj + fihja^a 

2 

+ [.9(''j>7«^ + .9*(''j>crj 



(9) 



where the counter-rotating terms have been ignored. The 
RWA is valid if the cavity field is nearly resonant with 
the atomic transition frequency, w « wq, and the coupling 
constant g{rj) is much smaller than lo, i.e., g{rj) <ti w. 

The master equation ^ can be solved numerically for 
various values of L and d and for different initial condi- 
tions. For our purposes here, it is sufficient to focus on 
the initial condition of the two qubits prepared in the 
spin anti-correlated state \s). We calculate the concur- 
rence between the atoms by tracing the density opera- 
tor Ps over the cavity field assumed to be in the vacuum 
state |0). 




FIG. 1: Concurrence as a function of the scaled time ujt and 
distance between the qubits evaluated without the RWA for 
go = u->, (uJ ~ uio)/ijJ = 0.01, k = O.Ilj and L — A/2. 



Consider first the evolution of an initial entanglement 
without the RWA. Figure [1] displays the time evolution of 
the concurrence for the initial spin anti-symmetric state 
|s) and for the case of a strong coupling of the qubits to 
the cavity mode, go = co. We see that the initial entangle- 
ment undergoes the discontinuity, the sudden death be- 
havior. The entanglement disappears quite rapidly over 
a very short evolution time, cot ^ I. The sudden death 
behavior appears for small distances between the qubits, 
where a strong coupling is present, g{rj) not much dif- 
ferent from go- For large distances, where g{rj) are very 
small, the discontinuity disappears and the entanglement 
decays asymptotically in time, the behavior predicted be- 
fore under the RWA. 

The entanglement sudden death seen in Fig. [T] seem 
puzzling at first, because it appears to contradict the 
predictions based upon a simple argument that with a 
single excitation present in the system, it is impossible 
to achieve the discontinuity as, according to Eq. ([3]), it re- 
quires population of the two-photon state. To resolve this 
problem, we now discuss the origin of the population of 
the two-photon state. We offer two complimentary views 
of the underlying physics. The first is provided by the 
derivation of the reduced density operator for the atoms 
using state vectors of the combined atoms plus the cavity 
field system. Assume that there is only a single excitation 
in the system that evolves under the RWA Hamiltonian. 
In this case, the space of the system is spanned by three- 
state vectors, |Ti)|0), |iT)|0) and where |0) and 
|1) are the zero-photon and one-photon Fock states the 
cavity mode. 

If we evaluate the reduced density operator p for the 
atoms by tracing over the cavity mode, we find 



P 



= Tr, 



P23(<)in>(Tii 



-P33WUT>(Til 

-P32{t)\m{iu 



(10) 
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where \i) refers to the cavity mode and pij are popula- 
tions {i = j) of the atomic states and coherences {i ^ j) 
between them. Thus, under the RWA, the evolution 
of the atoms does not involve the two-photon (upper) 
atomic state. 

Consider now the evolution of the system under 
the non-RWA Hamiltonian, that includes the counter- 
rotating terms aa~ and crja^. In this case, we must in- 
clude the processes that do not strictly conserve excita- 
tion number, when one of the atoms goes to the excited 
(ground) state by emitting (absorbing) a photon. Thus, 
the space of the system is now spanned by six-state vec- 
tors, |n)|0), |iT)|0), |n)|2), |iT)|2) and |TT)|1). 

This leads to the reduced density operator of the form 

P - Pii(i)Ui>ail+P22(i)|Ti)(iTI 

+ P33(i)UT>(Til+P44(t)|TT>(TTI 

+ P32(i)UT)(iT|+P23(i)|Ti>(Ti| 

+ Pi4(t)Ui)(TT|+P4i(t)|TT>(ii|. (11) 

It is evident from Eq. (fTT|) that the evolution of the 
atoms under the non-RWA Hamiltonian involves the two- 
photon state, that population of the state | It) becomes 
possible. 

A second explanation follows the principle of comple- 
mentarity between the evolution time At and uncertainty 
AE in the energy 

AtAE > H. (12) 

For the evolution time of the order l/w, that is considered 
here, the uncertainty in the energy is of the order hiu, the 
order required to achieve a non-zero population of the 
upper state | "ft) of the two-atom system. 

0.35 1 ^ ^ ^ ^ ^ 1 
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FIG. 2: Time evolution of the population of the two-photon 
state [4) calculated without the RWA for go = cj, (lu—ujo)/lj = 
0.01, K = Q.luj, L = A/2 and d = 0. 

To emphasize the participation of the two-photon state 
in the atomic dynamics under the non-RWA Hamilto- 
nian, we plot in Fig. [5] the time evolution of the pop- 
ulation of the two-photon state |4) = | It) for the same 



parameters as in Fig.[T]with d = 0. It is seen that initially 
unpopulated state |4) becomes populated in a compara- 
tively short time. 

As we have already pointed out, the two-photon state 
can be populated only if the counter-rotating terms are 
present in the interaction Hamiltonian. To show this ex- 
picitly, we look at the evolution of the entanglement un- 
der the RWA approximation. In Fig. [31 we illustrate the 
time evolution of the concurrence for the same parame- 
ters as in Fig. [U but under the RWA Hamiltonian 
We see that in contrast to the non-RWA case, the entan- 
glement evolves continuously without any discontinuity. 
Clearly, the counter-rotating terms are responsible for 
the entanglement sudden death. Under the strong cou- 
pling situation and small distances between the qubits, 
the entanglement oscillates with the Rabi frequency g^. 
The oscillations disappear and the evolution ceased for 
large distances, where the coupling constants are small. 
The fact that the RWA version of the interaction pre- 
cludes the entanglement sudden death suggests that cau- 
tion should be exercised in the studies of entanglement 
evolution in a strong coupling limit. We stress that the 
coupling strengths considered in this papers have not 
been realized yet. However, there are proposals of re- 
alistic systems involving hi gh Q cavities [2l|, [l^, or 
nanomechanical resonators [24], where such strong cou- 
plings could be realized experimentally with the current 
technology. 




FIG. 3: Concurrence as a function of the scaled time uot and 
distance between the qubits calculated with the RWA for go — 
LO, {uj — ujo)/uj = 0.01, K = O.ltJ and L = A/2. 

In summary, we have considered the time evolution 
of an entanglement initially encoded in a spin anti- 
correlated state of two identical qubits. What is special 
with this entanglement, that no discontinuity behavior 
has been predicted under the RWA. We have demon- 
strated the failure of the RWA in the description of the 
evolution of entanglement in the strong coupling regime 
of the qubits to the field they interact with. The results 
show that the better the RWA is in describing the dynam- 
ics, the weaker the entanglement sudden death behavior. 
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and only in the limit of a weak coupling, the RWA is an 
excellent approximation for the entanglement evolution 
of the spin anti-correlated states. 
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